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ABSTRACT 

We present a Herschel far-infrared study towards the rich massive star- 
forming complex G305, utilising PACS 70, 160 pm and SPIRE 250, 350, and 
500 pm observations from the Hi-GAL survey of the Galactic plane. The focus of 
this study is to identify the embedded massive star- forming population within G305, 
by combining far-infrared data with radio continuum, H2O maser, methanol maser, 
MIPS, and Red MSX Source survey data available from previous studies. By applying 
a frequcntist technique we are able to identify a sample of the most likely associations 
within our multi-wavelength dataset, that can then be identified from the derived 
properties obtained from fitted spectral energy distributions (SEDs). By SED mod- 
elling using both a simple modified blackbody and fitting to a comprehensive grid 
of model SEDs, some 16 candidate associations are identified as embedded massive 
star-forming regions. We derive a two-selection colour criterion from this sample of 
log (F70/F500) ^ 1 and log (F160/F350) ^ 1.6 to identify an additional 31 embedded 
massive star candidates with no associated star-formation tracers. Using this result we 
can build a picture of the present day star-formation of the complex, and by extrapo- 
lating an initial mass function, suggest a current population of « 2 x 10 4 young stellar 
objects (YSOs) present, corresponding to a star formation rate (SFR) of 0.01-0.02 
Mq yr _1 . Comparing this resolved star formation rate, to extragalactic star formation 
rate tracers (based on the Kcnnicutt-Schmidt relation), we find the star formation ac- 
tivity is underestimated by a factor of ^ 2 in comparison to the SFR derived from the 
YSO population. 
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1 INTRODUCTION 

The G305 star-forming complex is one of the most mas- 
sive and luminous star-forming regions in the Galaxy 
I Clark fc Ported . 12004 iHindson et all . |2010| ) , centred on the 
two optically visible open clusters Danks 1 & 2 (jPanks et all 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
f E-mail: a.faimali@herts.ac.uk 



1 1984 ). and the Wolf-Rayet star WR 48a. Located inside the 
Scutum-Crux arm within the Galactic plane at / = 305° , 
b = 0° , and found at a distance of ~ 4 kpc, it has a pro- 
jected dia meter of ~ 30 pc, and an estimated age of some 
3-5 Myr (|Clark fc Ported . l2004h . Coincident with G305 are 
numerous signposts of massive star formation, some of which 
are located on the periphery of the central cavity, such 
as infrared hotspots, compact and ultra-co mpact (UC) HII 
regions, H2O, OH, and m ethanol masers (|Urquhart et all 
120071 : IHindson et ail, l201Ch . 
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Figure 1. Three-colour (blue = 70 /mi, green = 160 /jm, red = 350 /mi) Herschel Hi-GAL image of the giant HII region G305. The 
70 |4m emission corresponds to hot dust emission 40 K), while the 350 um emission emanates from more cold dust 20 K). The 
positions of both Danks 1 & 2 (circles), and WR 48a (cross) are ovcrplottcd. 



We define a massive star as a star with sufficient mass 
to produce detectable radio emission associated with a HII 
region (i.e. M > 8 Mo). They play a key role in the Universe; 
their presence has a profound effect on the stellar and plane- 
tary formation process, while also on the physical, chemica l, 
and morphological structure of galaxies (|Kennicuttl 120051 ) . 
Their UV radiation output ionises the surrounding interstel- 
lar neighbourhood providing the principle source of heating 
in the interstellar medium (ISM), while also having an im- 
pact on subsequent star formation. Through various pro- 
cess, such as stellar winds and supernovae, their mechani- 
cal energy output serves as the energy and momentum in- 
puts to the surrounding ISM, helping to sculpt the struc- 
ture and energetics of the ISM, and hence in host galaxies 
|Zinnecker fc Yorkd. l2007l ) . 

Though they are important, in the shaping and evo- 
lution of their host galaxies, the physics of the formation 
and evolution of massive stars is unclear. Firstly, there is 
the issue of high dust extinction making direct observation 
difficult; sites of massive star formation tend to be embed- 
ded within very opaque cloud cores, which su ffer from vi- 
sual extinction greater than 100 magnitudes (|Ballv et all 
120051 ). Related to this also, is the fact that massive stars are 
predominantly foun d to form within dense, stellar clusters 
jde Wit et all 120051 ). In all, this makes for a particularly dif- 
ficult challenge in being able to study the origins of massive 
stars. 

Fig'tH underlines the dynamic morphology of the re- 



gion, with the central location Danks 1 & 2 and WR 48a 
thought to be the main source of energy input and the driv- 
ing force behind the expansion and clearing of the centra l 
diffuse HII region in the complex j Clark fc Ported . |2004| ). 
The suggestion is an interaction between the centrally em- 
bedded sources and the surrounding cloud, with an occur- 
rence of ongoing massive star formation being located in the 
hot dust emission sites (seen in blue in Fig[l} on the periph- 
ery of the central cavity. With the presence of numerous and 
different epochs of star formation in one location, and the 
relative close proximity, the G305 star-forming complex af- 
fords us an exceptional opportunity to study the nature of 
massive star formation, and investigate the environmental 
impact t his may have on the formation of future gene rations 
of stars lElmegreenl . 120021 ; lElmegreen fc Ladal , 1 19771 ) . 

This paper follows in a series of previous studies of 
massive star for mation within the G3 05 star-forming com- 
plex. Previously, iHindson et al.l (|2010l ) have focused on the 
reservoir for star formation within the region through ob- 
servations of NH3 emission, while tracing sites of active 
star formation through H 2O maser emission. To complement 
this, iDavies et alj (|2012| ) have recently carried out a near- 
infrared study of the two central open clusters Danks 1 and 
2; while current work has focused on identifying compact 
radio emission that is indic ative of UC HII regions towards 
G305 (|Hindson et al.l . l2012l ). 

In this paper, we present a far-i nfrared (far-IR ) stud y 
of the G305 complex using Herschel jPilbratt et all |2010| ). 
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in conjunction with radio continuum, H2O maser, methanol 
maser, MIPS, and Red MSX Source survey data, with the 
aim of identifying sites of embedded massive star formation. 
In this study we are able to resolve the embedded population 
within G305, since the Herschel far-IR observations are un- 
affected by dust, and combined with mid-IR data, constrain 
the luminosities of individual YSOs. By incorporating the 
luminosities of the embedded massive star-forming popula- 
tion with the initial mass function (IMF), we are able to 
determine the SFR of G305 and investigate the star for- 
mation history of the region. This resolved, Galactic SFR 
can then be compared to extragalactic SFR indicators to 
test whether the two regimes are consistent with one an- 
other, and identify where fundamental difference s may lie 
|Heiderman et all . l20ld ; lLada. Lombardi fc Alvesl l2010h . 

The present study serves as an example of how Her- 
schel data can be applied to Galactic star-forming regions, 
such as G305, in order to identify the high-mass stellar con- 
tent of such complexes, and how the star formation activity 
can be inferred from this population. Following from this, a 
comprehe nsive YSO counting approach, similar to that con- 
ducted bv lPovich et al.l (|201lf ) for the Carina complex, will 
be conducted to tackle the incompletenesses present at the 
intermediate, and low-mass range in a follo wing paper by the 
auth ors. By combining Sp itzer GLIMPSE jBeniamin et all 
120031 ). and VISTA VVV (|Minniti et all |2010| ) data along- 
side our present datasets we aim to conduct a complete cen- 
sus of the YSO population of G305, studying the physical 
parameters of these sources, their evolutionary stage, and 
spatial distribution within the complex. 

This paper is organised into the following sections; in 
Section 2 we present and discuss the details of the observa- 
tions and the data analysis, in Section 3 we present methods 
to identify sources detected in the far-IR through consider- 
ation of their spectral energy distribution, in Section 4 we 
present our discussion and place G305 in context with other 
star-forming regions, and in Section 5 we summarise our re- 
sults. 



2 OBSERVATIONS & DATA ANALYSIS 
2.1 Herschel Hi-GAL 

The Herschel Infrared GALactic plane survey (Hi-GAL) 
( Moli nari et al. I, l2010bl ) is an Open Time Key Project on 
board the ESA Herschel Space Observatory jPilbratt et all 
l2010l ), mapping a two degree wide strip of the inner Galac- 
tic plane, in the longitudinal range | / 1 ^ 60° and lati- 
tude ra nge |b| ^ 1°. The su rvey takes adv antage of the 
PAC S (|Poglitsch et al.l . l20ld ) and SPIRE (|Griffin et all . 
120101 1 cameras operating in parallel mode, scanning the sky 
in a raster fashion at a rate of 60 "s -1 , while PACS and 
SPIRE acquire data simultaneously. The survey has a to- 
tal angular resolution, 9, of 10" ^ 9 ^ 30" across five 
photometric bands at 70 and 160 /im with PACS and 250, 
350, and 500 /im with SPIRE. The overall aim is to cata- 
logue star-forming regions and study cold structures across 
the ISM. Using the broad spectral coverage available, the in- 
tention is to study the early phases of star formation, with 
particular focus on providing an evolutionary sequence for 
the formation of massive stars within the Galactic plane. 



The catalogue of compact infrared sources for the G305 
region is obtained using the CuTEx (CUrvature Th r eshold - 
ing Extractor) code highlighted by iMolinari et al.1 (|201lT ). 
The detection technique considers the curvature proper- 
ties of astronomical images, rather than source detection 
through signal intensity, by building a "curvature" image 
from the observed image using double-differentiation in four 
separate directions. The advantage of this approach is that 
resolved, compact sources are easily detected, while the dif- 
fuse thermal emission from cold dust associated with the 
fore/background is greatly reduced. Photometry estimates 
of candidate sources are then performed by fitting an ellip- 
tical two-dimensional Gaussian with an underlying planar 
inclined plateau. 

The completeness of the infrared source catalogue is 
estimated using a simulated field of some 1000 artificial 
sources, initially at 250 urn, comprising both a compact dust 
component (such as YSO envelopes, or dense molecular cores 
or clumps), and the diffuse emission present towards the 
Galactic plane jMolinari et al.l . [2010a). In all, source re- 
covery rate is found to be 90% for sources of peak fluxes 
of 0.2 Jy/pixel (equivalent to approximately 5<r), with peak 
fluxes being within 30 % of their input valu e in 80% of re- 
covered input sources <|Molinari et al.l . 120111 ). 

Identifying each source within the five different bands is 
obtained considering their basic positional association, start- 
ing with source extraction at the longest of the observed 
wavelengths, i.e. at 500 vim. From this first band, an asso- 
ciation is established with the next band (i.e. 350 vim) if a 
source is present within a search radius that corresponds to 
the Herschel beam size at the longer of the two considered 
wavelength, in this case a radius of 18.1" at 500 u.m. In the 
case where multiple associations are fou nd, the closest on e 
is kept. This is more frequent at 70 vim (|Elia et all |2010| ). 
At this stage the flux contribution at longer wavelengths has 
not been split up into contributions from multiple associa- 
tions at 70 urn, however we find 87% of 500 vim sources are 
associated with a single 70 vim detection. Of the remaining 
500 vim sources with multiple 70 vim associations, the major- 
ity are found to have the bulk of the 70 vim flux assigned to 
their primary 70 vim counterpart. We find only 3% of 500 vim 
sources to have lost a fraction of 70 vim flux through multi- 
ple associations. This loss is deemed negligible, being on the 
order of 1.5% of the total 70 vim flux. 

Finally, the celestial coordinates assigned to the sources 
are those which correspond to the shortest wavelength as- 
sociation, which by defini t ion wi ll have the highest spatial 
resolution C Molinari et al.l . boiOal ). As a result of this, in to- 
tal some 3288 infrared sources are detected across the G305 
region in Hi-GAL; of these some 1913 sources are detected 
at 70 vim, 1658 at 160 vim, 1257 at 250 vim, 856 at 350 vim, 
and 530 sources at 500 vim. Not all sources are detected in 
all wavebands, and this can be accounted for by either po- 
sitional association failing, or that the flux obtained for a 
source in a particular band had been corrupted by factors 
such as source crowdin g, or confused background conditions 
|Molinari et all 1201031 ). 

2.2 Ancillary Data 

To complement our data set we utilise mid-IR, radio, H2O 
and methanol maser data matched to the Hi-GAL far-IR 
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Figure 2. Distribution of P-statistics of identified associations 
(solid blue line), against the distribution of P '-statistics obtained 
from Monte Carlo associations (dashed red line). 



and UC HII regions lo cated throughout the Galactic plane 
(|Urquhart et all 120081 ), and the 24 um point source cata- 
logue f rom the Spitze r MIPS Galactic Plane Su rvey (MIPS- 
GAL) (|Mizuno et all 120081 ; ICarev et all . l2009f ) which aims 
to identify all high-mass protostars located in the inner 
Galactic disk. 



2.3 Identification Procedure 

Our aim is to identify a statistically reliable sample of coun- 
terparts to the Hi-GAL infrared catalogue, by matching this 
catalogue to the ATCA radio, MMB, RMS, Mopra H 2 
maser, and MIPSGAL ca talogues. By usin g a frequentist 
technique, highlighted by I Lilly et al.l (|l999l ). we can iden- 
tify all associations found within G305, by establishing the 
probability that matches to Hi-GAL sources are not the re- 
sult of a chance alignment. We approach this using a Monte 
Carlo method, with the technique outlined in the following 
subsection for the case of Hi-GAL infrared and ATCA radio 
counterparts. 



data to broaden our view of the G305 complex, aiding in the 
consideration of the morphology, and star formation within 
the region. 

We firstly make use of the 5.5 GHz radio continuum 
observations towards G3 05 using the Austral ia Telescope 
Compact Array (ATCA) (|fflndson et al.l , l2012h as a means 
to distinguish UC HII regions across the complex. A key 
method of detecting HII regions is through the associated 
radio emission; at wavelengths greater than 3 mm, thermal 
dust emission drops off, and the dominant detection is due 
to thermal free-free radio emission from the ionised gas in 
the inner region surrounding the massive star. The Lyman 
continuum photons that are emitted from the central star 
are balanced by recombination within the volume of ionised 
hydrogen, producing the detected radio emission. For this 
emission, stellar models suggest stars of spectral type B3 
or earlier are capable of producing the required ionising 
flux of Lyman cont i nuum photons (E>13.6eV, i.e. 912 A) 
(|Crowther fc Contll I2003T ) . Such a spectral type is selected 
since the Lyman continuum flux then begin s to rapidly dro p 
off with decreasing effective temperatures (|Panagial . I1973T ). 

To also locate areas of ongoing star formation, we em- 
ploy both 22 GHz H2O maser observations using the Aus- 
tralia Telescope Nati onal Facility (ATNF) telescope Mopra 
(|Hindson et al.l,l2010l ). and the 6.7 GHz methanol maser ob- 
servations of the Me t hanol Multi-Beam survey l| Green et al.l . 
120091 ; I Caswell et all, I2010T ) with the Parkes 64 m radio tele- 
scope and follow ups by the ATCA. Since their discovery, 
methanol masers have been recognised as one of the mos t 
distinct signposts of massive star formation (lMentenl.ll99lf) . 
and o nly found close to high- mass young stars ( Minier et al.l . 
120031 ). H2O maser emission, on the other hand, has also been 
shown to provide a useful indicator of both low- and high- 
mass star formation, and are found to associ ate mainly wit h 
both hot molecular cores and UC HII regions (jFuruvalEiblih . 

Finally, we also make use of mid-IR data from both the 
Red MSX Source (RMS) survey that aims to identify a large 
sample of genuine massive young stellar objects (MYSO) 



2.3.1 Association Probabilities 

Starting from the matching of both the Hi-GAL infrared 
and ATCA radio catalogues, we need to consider the pos- 
sibility of chance alignments and provide an estimate as to 
how reliable each individual match may or may not be. We 
therefore need to derive some statistical argument that con- 
siders the probability that a candidate compact radio source 
is indeed correctly identified within the search radius of the 
associated infrared source. A method b ased on the pos i tiona l 
coin cidence, similar to that adopte d bv lDownes et all (|1986T) 
and ISutherland fc Saunders! (1 19921 ) is employed. 

The necessity to work out a statistic which can inform 
whether an association is indeed true, or by chance can 
be shown by considering the ATCA radio data. A certain 
number of the 71 radio sources identified across the field 
will be accounted for by extragalactic background sources. 
We can estimate the total amount of background sources 
e mpirically from the ex tragalactic source count approach 
of lAnglada et all (|l998h . The number of background radio 
sources, ( N Radio ) , that would be expected to be observed 
given the region size, and frequency of observation is defined 
by: 



|q i.is -°- 7! 



(i) 



where 9f is the diameter of our observed field, 9a is the 
FWHM of the primary beam (in arcminutes), and So is the 
sensitivity of the radio observations at 5.5 GHz. Using this 
approach, we find that some 60 ± 8 backgrou nd sources in to- 
tal sh ould be detected across the G305 field (|Hindson et al.l . 
120121) . 

As an example, to explain the positional coincidence 
method, we use our identified infrared-radio associations. 
For each association we calculate that for each ATCA com- 
pact radio source candidate with magnitude m at a distance 
r from the matched to the Hi-GAL infrared source, there is 
a surface density N a of radio sources brighter than m across 
the G305 field. We therefore obtain the mean number ft of 
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Table 1. Optimum matching radius for each data set to the Hi- 
GAL G305 field. 



Data Set 


True ID 


False ID 


Optimum Matching Radius ( ) 


ATCA 


11 


2 


IS 


MMB 


14 


2 


15 


RMS 


9 


> 


5 


H 2 Masers 


4 


1 


15 


MIPS 


500 


63 


10 



chance sources that are closer and brighter than the matched 
to candidate ATCA radio source: 



/' : 



(2) 



We utilise this statistic to suggest the fraction of sources 
within a sample size of n, that we would expect to have an 
incorrect candidate association identified within a matching 
radius r. In a sense, the probability of the association is: 



P — [1 — exp(— /x)] « jj, for (i<l 



(3) 



From this, we find that if the P-statistic, P^Cl for an indi- 
vidual source, then that particular identified association is 
unlikely to be the result of a chance association. 

However this statistic alone does not provide us with a 
firm argument to the reliability of associations; merely the 
chance that the particular source would have an association 
within the specific matching radius. Rather, a more reliable 
means to measure whether an individual association can be 
deemed correct, is to compare the number of identifications 
in the total sample with a particular P-statistic, against the 
amount of associations, nP, that would have been expected 
from a randomised association between the infrared and ra- 
dio populations. 

To do this, we run 5000 Monte Carlo simulations that 
follow the previous matching criteria between both Hi-GAL 
and ATCA catalogues; where the positions of each com- 
pact radio source have been randomized, with the only con- 
straint being that they remain within the field. From these 
Monte Carlo results, we can then calculate the P '-statistic 
for each of the identified Monte Carlo candidate matches. 
By comparing the P-statistic of our associations to those 
P '-statistics from the randomised sample, we examine the 
ratio of sources with a particular value of P to the similar 
value of P' from the spurious identifications. If this ratio is 
high, i.e. there are many associations with a particular value 
of P compared to that of P' , then we can mark that partic- 
ular association as a secure identification. This is shown for 
a matching radius of 30" in Fig. [2] 

2.3.2 Separation Cut-off 

A useful outcome from the use of P-statistics is that it al- 
lows for the determination of the separation cut-off for each 
dataset; essentially the maximum matching radius at which 
one considers the majority of associations found to be reli- 
able. The selection of the separation cut-off depends on sev- 
eral factors. Clearly a smaller cut-off gives rise to fewer coun- 
terparts, increasing the likelihood of missing a true counter- 
part. Conversely, a larger cut-off increases the risk of associ- 
ating to a very bright but unrelated counterpart. From the 
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Matching radius (arcseconds) 

Figure 3. Distribution of radial offsets between Hi-GAL sources 
and ATCA radio counterparts (open histogram). The dashed line 
gives the expected cumulative number of false radio IDs as a 
function of the separation cut-off, while the dot-dashed line shows 
the cumulative number of true radio IDs for each separation cut- 
off. In this context, the optimum matching radius is found to be 
15". 



P-statistics we are then able to derive the optimum match- 
ing radius for th e data using the technique highlighted by 
iDve et all (|2009h . Firstly, we match between the Hi-GAL 
and ATCA data at a minimum matching radius of 1". By 
then comparing the P-statistic of our associations to those 
P'-statistics from the randomised sample, we examine the 
ratio of sources with a particular value of P to the similar 
value of P' from the spurious identifications. If this ratio is 
high, i.e. there are many associations with a particular value 
of P compared to that of P' , then we can suggest that par- 
ticular association as a secure identification. This process is 
stepped through, at intervals of 1" up to a maximum match- 
ing radius of 30" . As a result of this we are then able to 
compute a true and false ID rate as a function of separation. 

The result is shown in Fig. with the true ID rate, 
the false ID rate, and the separation between Hi-GAL and 
ATCA true associations in the G305 field shown as a func- 
tion of the separation cut-off. From this we can see that 
the true ID rate levels out at 15" (with a total of eleven 
true associations), and only increases again at 25" onward 
(with a further three true associations being found). These 
additional three true associations found at radii ^ 25" cor- 
respond to Hi-GAL matches to the brightest ATCA radio 
sources found in the G305 field. In choosing the search ra- 
dius, we wish to maximize the number of secure, unambigu- 
ous identifications and to minimize the number of real coun- 
terparts missed. It is not seemingly clear as to why there is 
a well-defined P-statistic peak at a lower matching radii; 
being due to faint, real radio identifications, contesting (in 
terms of low P-statistic) with rare, brighter sources, which 
are unrelated to the Hi-GAL source (such as background 
radio galaxies). By increasing our matching radius, any cor- 
rect identifications will tend to drop off, with the levelling 
out in true identifications being balanced by those bright 
unrelated radio sources present; since they are very bright, 
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their P-statistic tends to zero as soon as they are within the 
matching radius. Therefore allowing a matching radius well 
beyond 15" is not suggested, as the amount of secure associ- 
ations tend to level out, while matches to contaminant radio 
sourc es will preserve the number of perceived true associa- 
tions |lvison et all . 120101 ). 

We also note from Fig. [3] that the false ID rate only 
becomes noticeable for 15" onwards, with two possible false 
IDs found at 15". Finally, the separation between Hi-GAL 
and ATCA true associations, shown in the open histogram, 
emphasises the fact that the majority of true associations are 
found within a matching radius of ^ 15". The outcome of the 
matching procedure returns an optimal separation cut-off of 
15" for the radio counterparts in the Hi-GAL field. 

Suggesting that a radio source some 15" from an IR 
source is a secure match may seem counter-intuitive, but 
two factors need to be taken into consideration; the Her- 
schel beam size, and the physical nature of the IR sources. 
Our goal is the identification of a sample of embedded mas- 
sive star-forming regions that themselves are extended in 
nature, on average aroun d 29" in size, and with the beam 
size some 18" at 250 urn (|Traficante et all 1201 if ) . a secure 
match within 15" is acceptable. 

The approach is repeated for the MMB, RMS, H 2 
Masers, and MIPS data sets, with the total number of as- 
sociations found and the optimum matching radius shown 
in Table [T] The result of the technique is a statistically ro- 
bust sample of true associations that have been identified 
over multiple datasets, which can later be investigated and 
analysed. 



3 SPECTRAL ENERGY DISTRIBUTIONS OF 
SOURCES 

To obtain estimates of the physical properties of our asso- 
ciations, the observed SEDs were fitted with firstly a sim- 
ple modified blackbody, incorporating the Herschel Hi-GAL 
observations at 70, 160, 250, 350, and 500 urn. We then can 
obtain estimates of the bolometric luminosity of each as- 
sociation, by broadening the SED coverage with both the 
MSX 21 and MIPS 24 urn fluxes, and fitti ng to the grid 
of SED models from iRobitaille et all | |2006j L We find that 
the modified blackbody fails to reproduce the observed flux 
at A 70 vim, hence associations with solely Herschel de- 
tections are best suited to modified blackbody fits. Those 
associations with detections at 21 and 24 um, suggesting a 
warm embedded YSO, are more reliably reproduced with the 
IRobitaille et al.l l|2006r l SED models, that incorporate both 
the central embedded source, and the surrounding envelope. 

3.1 Modified Blackbody Fitting 

In order to firstly be able to derive the basic physical prop- 
erties for each of our associations such as dust temperature 
and emissivity index, it was necessary to model the SEDs 
with a modified blackbody. We are justified in taking the ap- 
proach of a simple modified blackbody, since the wavelength 
coverage measured by both PACS and SPIRE trace the peak 
of the dust SED. Far-IR emission is due to large dust grains 
(15-100/im), which are more stable and tend to dominate 
the total dust mass, while also tending to trace all phases of 
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Figure 4. Fits to two associations using a modified blackbody 
SED. Physical properties derived from the fit are shown, along 
with minimum \ 2 for the best fit. The red solid line represents the 
best fit model, while the grey solid lines represent other models 
with a good fit to the data. 



the ISM. Therefore, measurements in Hi-GAL will be most 
sensitive to temperature variations, while also providing an 
accurate tracer of the overall ISM column density. 

Since we have a embedded massive star within a dust 
cocoon we cannot assume a blackbody, but rather a modified 
blackbody that takes into consideration the optical depth, 
and the dust emissivity, as clearly emission is not optically 
thick at all frequencies. Thus w e adopt a modified blackbody 
fit to the SED of the form (|Ward-Thompson fc RobsorJ . 
Il990h : 

F v = f2i?„(T)(l — e~ T ) [Jy] (4) 

where f2 is the effective solid angle of the source, in sr, B„(T) 
is the Planck function, T is the dust temperature, and r is 
the optical depth such that: 
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(5) 



where j3 is the dust emissivity index, and v a is the critical 
frequency at which r = 1. 

The free parameters that are derived from the fit in- 
clude the dust temperature, the dust emissivity index, and 
the critical frequency. We also note that throughout the cal- 
culations each Herschel flux was assigned an uncertainty of 
20%, this being the calibration error. 

Fitting was performed via a \ 2 minimisation, by consid- 
ering the observed flux at each of the five Hi-GAL wavebands 
available for every i ndividual associatio n. The \ 2 minimisa- 
tion was such that (|fiunter et all I2OO0I ') : 



x 2 = E 



F v n , model 
.observed 



(6) 



where the ratio of the model flux to the observed flux is used 
to give equal weighting to each different wavelength regime. 

With these parameters we can then derive the total 
(gas + dust) mass for each association using both the Her- 
schel fluxes at 500 /im, and the corresponding free param- 
eters derived from the minimum x 2 fit- We note that dust 
masses were calculated with an opacity rsoo/jm , since at this 
wavelength all sources were found to be optically thin, al- 
lowing for us to sample the dust at all depths. If dust masses 
were calculated at shorter wavelengths, where some sources 
were found to be optically thick, we would simply be trac- 
ing the dust distribution in the outer layers, and not in fact 
determining a total dust mass value. Therefo re, by following 
the method highlighted in iHildebrandl l| 19831 ): 



f n 2 



where D is the source distance, and C v is the mass coefficient 
(a factor that combi nes dust opacity and gas-to-dust ratio) 
dKerton et all l200lt ): 



C v = 



MaK v 



(8) 



where M g is the gas mass, Md is the dust mass, and k v is 
the dust opacity. A value of the mass coefficient of C„ = 
50 gcm~ 2 at 850 um w as initially chosen from those quoted 
in lKerton et al.l (|200ll ). assuming a gas to dust ratio of 100 
and a dust emissivity index /3 of 2. This value was then 
scaled to a value at 500 /im, by adopting the value of j3 from 
the minimum \ 2 fits, and assuming that: 



C v oc v 



(9) 



Fig.|4]shows the result of the modified blackbody model 
SED fit for a sample of associations, with both the best fit 
model and thos e models deemed go od fits from the criterion 
X 2 I Ndata ^ 2 jPovich et all l201ll ) shown. Included are the 
corresponding free parameters of the minimum \ 2 fit, being 
the dust temperature, emissivity index, and source opacity 
at 500 ^m (rsoo^m)- 




Figur e 5. Fits for two associations using the iRobitaille et al] 
1120071 ) fitting technique. Black solid line represents the best fit 
model, with the grey solid lines showing all other models pro- 
viding a good fit to the data. The black filled circles show the 
Hi-GAL and MSX/MIPS data with error bars shown. 



3.2 SED Fitting 

We next obtain estimates of the bolometric luminosity by 
fitting the observed SEDs w ith the grid of young stellar 
object (YSO) model SEDs of IRobitaille et al.l (|200rj ). Since 
the SEDs of embedded YSOs tend to pea k at 100 urn ,we 
requir e photometric data at A ^ 10 \im ijMottram et all 
l2011bl : |Pwich et all [ioTlh . hence to broaden the SED cov- 
erage we incorporate the 21 and 24 [im fluxes obtained from 
MSX and MIPS. For these sources, the standard modi- 
fied blackbody adopted earlier will fail to produce accurate 
SEDs at wavelengths A ^ 7 ^m, thus a mu lt i-com ponent fit 
using the model SEDs of IRobitaille et al.l (l2006t) that are 
then fi t with the online SED fitting tool of IRobitaille et all 
(l2007h. that ar e based on the YSO/disk/envelope models of 
IWhitnev et all (|2003h a re employed. 

The model grid of IRobitaille et al. I (|2006l ) consists of 
some 200,000 model SEDs incorporating a vast range of 
possible evolutionary stages, from an embedded protostellar 
phase to pre-main-sequence stars with low-mass circumstel- 
lar disks. Fitting of models to the data is done by varying 
the visual extinction, Ay, for a number of distances d be- 
tween dram - d max , in order to determine the optimum SED 
model and parameters set. A range of 0-20 is selected for 
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the visual extinction (|Leistra et al.1 . 120051 ). while an averaged 
distance of 3.2 - 4.4 kpc derived from both the kinematic dis- 
tance to G305 and the spectr oscopic distance of Danks 1 & 
2 is used (|Davies et al.ll201.2f ). The optimum SED model is 
determined through linear regression, with the result shown 
in Fig. [5] where all models that fitted with a \ 2 value satis- 
fying x 2 ~ Xbest < 3 x n data shown also. 

The majority of sources, with no 21, 24, or 70 urn emis- 
sion, have fitted SEDs that peak at wavelengths A 160 um 
and are found t o have an averaged temperature of w 14 K. 
iBontemps et al.l (|2010T ) characterise YSOs within the Aquila 
rift complex from pre-stellar sources on the presence of ei- 
ther a 24 or 70 urn counterpart, with this emission originat- 
ing from warm dust within the inner regions of the YSO 
envelope. Moreover, such emission is unlikely the result of 
external heating from the interstellar radiation field produc - 
ing a detectable 70 uxn counterpart (|Giannini et al.l . l2012f ) , 
making sources with 21, 24, or 70 urn emission likely pro- 
tostellar sources. Indeed, iDunham et al.l l|2008l ) show that 
70 urn emission is closely correlated to the overall luminos- 
ity of protostellar sources. Using this selection criteria, the 
suggestion is that these identified sources correspond to an 
earlier stage (possibly pre-stellar) that are poorly fitted by a 
protostar embedded in a dust envelope, and are more suited 
to a simple modified blackbody fit, yielding values of dust 
mass, temperature, and emissivity index. 



4 DISCUSSION 

In this section we discuss the global distribution of the em- 
bedded massive star-forming population within the G305 
complex, and the general method of identification of this 
population; through the combination of the SED morphol- 
ogy and derived parameters. Using these identified sources, 
we are able to comment on the present-day star formation 
activity of G305, and place it into the context of Galactic 
star formation. 



4.1 General Properties Of Sources Within G305 

We firstly consider the global properties of G305 ob- 
tained from SE D fittin g, by both modified blackbody and 
iRobitaille et al.l (|2006h YSO models, for each association 
identified within the G305 region. For our fitting we con- 
sider sources with detections of Ndata ^3 as acceptable to 
fit an SED to, since we have three free parameters in our 
modified blackbody fitting, sources with limited detections 
will unlikely produce a reliable SED; this leaves us with 
some 503 sources. From these fits , we deem sources with 



X 2 /N d 



,2 (|Povich et al.l . 1201 J ), as having a reliably fit- 



ted SED, yielding a total of 359 well fitted sources. 

In Fig. [H] we present the distribution of the modified 
blackbody parameters of the dust temperature, dust emis- 
sivity index, and the dust temperat ure, along with the bolo - 
metric luminosities obtained from IRobitaille et al.l (|2006T l 
model SED fitting. We find that dust temperatures lie within 
the range of 10 to 42 K, with a median value of ~ 21 K. 
For the dust emissivity index, a range between 0.8 to 2.8 is 
found, with the median of f» 1.8. The bolometric luminosity 
is sampled between 10 to 10 4 Lq , with a median of » 300 Lq , 




Figure 6. Distribution of parameters for all sources with a re- 
liable SED fit within G305 (solid black line), compared to the 
16 identified candidate embedded massive star-forming regions 
(red dot-dashed line). Top left: dust temperature. Top right: 
dust emissivity index. Bottom left: bolometric luminosity. Bot- 
tom right: dust mass. 



while the dust mass lies between 1 to 10 4 M 
of ^540 M . 



with a median 



4.2 A Far-IR Selection Criterion For Embedded 
Massive Star Formation 

Previous studies that have attempted to distinguish a sam- 
ple of MYSOs ha ve relied on colour selection criteria based 
on IRAS d a ta dPalla et al.l. Il99ll; [Molinari etafl, Il996l ; 
IWalsh et all . Il997l ; ISridharan et all 120021 ). however these 
are restricted by selection effects, limitations in sensitivity 
and resolution, and also their incomplete Galactic coverage. 
The majority of these studies are biased in that they rely on 
bright IRAS sources, that with a resolution of 5' at 100 urn , 
are often found removed from the dense, confused regions 
along the Galactic plane where the majority of star forma- 
tion is expected to be found. The result of these studies is 
a restricted sample that, due to selection issues, may not 
entirely represent the general MYSO population, making 
any extension to study massive star- formation problematic 
jUrquhart et al.1, 12008T ). 

It has been suggested that both the SED morphology, 
and the bolometric luminosity can prove effective in deter- 
mining a sample of em bedded massive star-forming objects 
jMolinari et alll2008aD . With the bulk of the YSO emission 
being in the fa r -IR po rtion of the SED, due to cold dust 
|Molinari et all I2008H ) . the Herschel Hi-GAL observations 
are ideal in accurately constraining the luminosity. 

From the candidate associations identified, those 
sources found to have either a radio, MMB, H2O maser, or 
RMS counterpart are known sites of massive star formation; 
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Figure 7. Location of candidate embedded massive star-forming regions and their relative association overlaid onto a 160 fim greyscale 
Hi-GAL image; ATCA radio data (blue circles), MMB methanol masers (cyan asterisks), RMS sources (red crosses), H2O masers (purple 
diamonds). Numerous regions are found to be associated with multiple star formation tracers. 



tending to be luminous sources (i.e. > 1O 3 L.0). Using this 
sub-sample we are able to identify a sample of known em- 
bedded massive star-forming regions from the physical prop- 
erties derived from, and morphology of, their SEDs (shown 
to peak at 100 um, as shown in Fig. [5]). We refer to this 
population of Hi-GAL sources that are associated with ra- 
dio, MMB, H2O maser or RMS counterparts as the embed- 
ded population. Added to this, we employ a selection cut 
in the bolometric luminosity of 10 3 Lq , that corresponds to 
the minimum spectral type that we define as a massive star 
(i.e. M > 8M0). In total we find some 16 candidate embed- 
ded massive star-forming regions, that match these selection 
criteria; their respective properties and relative associations 
are shown in Table [2] We also note the location of these as- 
sociations, shown in Fig. with the majority being located 
in the hot dust locations along the periphery of the central 
cavity, suggesting an interaction between the central sources 
of Danks 1 & 2, and the surrounding material. 

Based on this sample of 16 embedded massive star- 
forming regions, we are able to derive a two-colour selec- 
tion criterion to identify the overall embedded population 
of the G305 complex. We find that the 70-500 \lm and 160- 
350 um colours are most sensitive to the embedded popu- 
lation, shown as asterisks in Fig. [S] As can be seen from 
Fig. [8] the embedded population are tightly confined to 
one area of the colour plot, and can be distinguished from 



the remaining G305 population, shown as circles. From this 
we suggest a far-IR colour selection criterion for embed- 
ded massive star- forming regions of log (F70/F500) ^ 1 and 
log (F160/F350) ^ 1.6, yielding a further 31 embedded mas- 
sive star candidates with no associated emission and lumi- 
nosities > 10 3 Lq, as shown in Fig. [8] From these 31 candi- 
dates we find the faintest source to have a peak 70 |J.m flux 
of 1.02 Jy/pixel, and compared to the 90% recovery rate of 
sources discussed in Sect. 2.1, we do not expect any further 
more deeply-embedded massive star-forming regions to be 
found within G305. 

Currently, the nature of these 31 candidate embedded 
massive star-forming regions is unclear; from Fig. [8] these 
candidates are predominately found with bolometric lumi- 
nosities of w 10 3 L0 however at least 3 sources are found a 
luminosities ^ 10 4 Lq with no corresponding tracer of mas- 
sive star formation. A lack of association to ATCA radio 
sources could be accounted for by localised noise found to- 
wards bright large-scale emission present with G305, confus- 
ing possible associa t ions to real compact emission present 
(see iHindson et al.l (|2012l ) Sect. 2.2 for a detailed discus- 
sion of the ATCA data reduction process). Another pos- 
sibility could be t he strong variabilit y of both methano l 
l|Green et al.l . l2012t ) and water masers (jBreen et al.l . l20ld ). 
to such an extent that they display no common features at 
the present epoch. Aside from the possible reasons for lack of 
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Table 2. Derived physical properties for all identified embed- 
ded ma ssive star- forming regio ns, from both modified blackbody 
fits and lRobitaille et aljT2007f> SED fitting technique, along with 
found associations. 



Hi-GAL 


P 


T 


T 500(im 


M Dusf 




Association(s) 


Source Index 




(K) 


(io- 3 ) 


(M ) 


(10 3 L ) 




938 


1.0 


35 


20.1 


4300 


48.6 


MMB, H 2 maser 


948 


1.1 


29 


15.4 


3900 


7.95 


RMS 


972 


1.2 


41 


11.4 


1200 


20.8 


ATCA, MMB 


1184 


1.1 


25 


27.1 


4700 


11.1 


MMB 


1800 


1.4 


35 


29.4 


3500 


60.8 


ATCA, MMB, RMS 


1804 


1.0 


32 


24.8 


5200 


17.4 


ATCA 


2114 


1.7 


24 


1.29 


200 


2.91 


RMS 


2153 


1.2 


29 


9.14 


3700 


4.13 


MMB, RMS, MIPS 


2212 


1.1 


23 


14.3 


4400 


4.68 


RMS, MIPS 


2363 


1.1 


27 


14.2 


3200 


12.2 


ATCA, RMS 


2383 


1.0 


36 


20.0 


1700 


57.7 


ATCA, MMB, RMS 


2627 


1.4 


28 


4.18 


1500 


3.65 


MMB 


2902 


1.0 


37 


42.2 


4200 


48.6 


MMB 


2923 


1.1 


24 


16.1 


3700 


3.24 


MMB 


2994 


1.0 


_><) 


23.1 


3900 


7.95 


MMB, MIPS 


3032 


1.4 


28 


4.57 


1500 


3.75 


MMB, RMS 



associated tracers, this sample of candidates may also sug- 
gest an earlier, very young embedded phase present within 
G305. The possible nature of these additional candidates is 
particularly interesting, and warrants further investigation 
at a later date. 



4.3 The Present-day Star Formation Rate Of The 
G305 Complex 

With a sample of embedded massive star-forming objects, 
one opportunity that is available is to study the recent star- 
forming history of the complex. Determining the star forma- 
tion rate (SFR) at a local level is crucial in determining the 
global Galactic SFR, helping to unveil any mechanism s that 
may lead to global scaling laws l|Molinari et all l2010b[ ) . The 
SFR, along with the IMF, express the population of mas- 
sive stars within the Galaxy, and determine what the impact 
on the local environment is, such as the composition of the 
ISM, the rate of feedback fro m massive stars , and t he rate 
of conversion of gas into stars (|Calzetti et all l2010al ) . 

Given our position within the Galactic disk, direct SFR 
indicators using optical/UV tracers will fail to reproduce 
an accurate SFR due to high extinction rates of the dusty 
ISM. However, far-IR observations, unaffected by extinction, 
provide us with the ability to resolve the YSO population 
associated with HII regions, allowing constraints on the IMF 
and stellar ages, yiel ding a detailed star forma tion rate of 
Galactic HII regions (|Chomiuk fe Povichl, [20lTI ) . 

A SFR derived from a resolved, YSO counting ap- 
proach, or from that inferred from the total infrared lu- 
minosity has to assume a 'steady-state' approximation to 
reliably trace the star formation activity of the region. The 
assumption in these calculations is that both the rate of 
massive star formation, and the rate that massive stars 
evolve off the main sequence , is in approximate equilibrium 
l|Krumholz. Dekel fc McKeel , I2012T ). For this to be true, the 
requirement is that the age of the region be longer than 
that of the UV e mitting population used to trace the SFR 
(|Kennicuttl Il99&t h We show below that this is the case for 
a realistic star-forming timescale. 
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Figure 8. Colour - colour plots of all Hi-GAL sources found 
in the G305 field (blue circles) , and known sites of massive star 
formation (red asterisks). Dashed lines indicates the boundary of 
the region used for distinguishing sites of embedded massive star 
formation from other sources, at a luminosity of > 10 3 Lq. 



4-3.1 The Embedded Massive Star Formation Rate 

With our sample of 16 identified embedded massive star- 
forming regions, and the further 31 embedded massive star 
candidates found, we are able to comment on the present 
star formation history of the G305 complex. If we make the 
assumption that for each embedded massive star-forming re- 
gion identified, the most massive star present produces the 
majority of the bolometric luminosity and is also accompa- 
nied by a cluster of lower mass stars, we are able to scale 
the IMF accordin gly; for our calcu lations, we adopt a simple 
IMF proposed bv lSalpeterl (|l955l) . By comparing the calcu- 
lated bolometric luminosity for each candidate region found, 
to the luminosi t ies calc ulated for a sample of MYSOs from 
iMottram et~aLl (|2011al ). we are then able to estimate the 
most massive star for each region. 

In order to extrapolate the IMF from observed mass 
to lower mass, we are obliged to select both a lower and 
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upper mass limit that all observed YSOs fall within the se- 
lected range. For our purposes we adopt a lower mass of 
0.1 Mq and an upper l imit of 50 Ma, as used when calculat- 
ing the Galactic SFR (|Robitaille fc Whitnevl l20ld ). 

By adopting a Salpete r IMF, which best s ampl es the 
high-mass tail of the IMF (jZinnecker fc Yorkd . l2007h . and 
assuming a constant SFR therefore, we arrive at « 10 4 
YSOs present, which corresponds to a total mass in stars 
of « 8 x 10 3 M . Since we consider the total mass in stars, 
we select a typical timescale for that mass to assemble, which 
would be the time to rea ch the pre main-sequence of 0.5 Myr 
(jOffner fc McKed . l201lf ). from this we attain a SFR for 
G305 of « 0.01 -0.02 M Q yr -1 . In this scenario, the 'steady- 
state' approximation s hould hold since the timescales of 
both MYSOs, « 1 4 vr dMottram et all 1201 left , and UC HII 
regions, w 10 5 yr jComeron fc Torral . 119961 ). are some 5-50 
times shorter than our assumed timescale for the star for- 
mation within G305 of ~ 0.5 Myr. 

A similar approach has been taken in iHindson et all 

(|2012h . using the identified UC HII population (some five 
identified in total) of G305 to derive a SFR, over the last 
0.5 Myr, of 0.002-0.004 M yr" 1 . This rate is considered a 
lower limit due to the incompletenesses in the ATCA radio 
data, with a uv cut implemented to emphasise the com- 
pact, small scale radio e mission assoc i ated with UC HII 
regions. For comparison, iDavies et all j|2012h . have deter- 
mined the SFR for the G305 complex over the last 5 Myr 
from the calculated age and mass of the two central open 
clusters Danks 1 & 2. The SFR was found to be « 0.002 - 
0.005 Mq yr" 1 , which is comparable to our derived SFR 
from counting the embedded YSO population. However, tak- 
ing our derived SFR of ^0.01-0.02 M Q yr -1 , it is clear 
that the star formation activity of G305 has not remained 
constant over the last 5 Myr. If this were to be the case, 
then we would expect to observe some 75 000 Mq of stars 
to be observed within the complex, which is entirely not 
the case. Th e 'collect and collap s e' mo del of star formation, 
proposed bv lElmegreen fc Ladal (| 19771 ). requires a period of 
time after the formation of the central ionising source(s) 
for material within the surrounding molecular clouds to be 
swept-up by the expansion of the HII region. If this expan- 
sion continues for a sufficient time, the swept-up shell of 
material becomes self-gravitating and is then expected to 
fragment, entering a phase of collapse possibly leading to 
th e formation of new stars. Followin g the approach taken 
bv lDale. Bonnell fc Whitwortbl |2007l ). who have calculated 
this 'fragmentation time-scale' for a uniform molecular cloud 
of pure molecular hydrogen, we take the ionis i ng ph oton flux 
of the G305 complex from IClark fc Porter! (|2004h to esti- 
mate a similar timescale. We find that after the formation 
of Danks 1 & 2, there would be a delay of ~ 2.4 Myr un- 
til the next generation of star formation occurred within 
the complex. This would support the scenario that the star 
formation within G305 was not continuous, but more likely 
characterised by punctuated star formation over the lifetime 
of the complex. 

We find also that our derived SFR for G305 is compa- 
rable to other well known massive star-formi ng complexes 
in th e Galaxy, n amely the Carina comple x (|Povich et all 
|2011| ). and M17 (|Chomiuk fc Povicbl. [201ll ) . We stress that 
the derived SFR value is based on a small sample of high 
mass stars, and has been extrapolated over a large range 



of stellar masses; when considering the lower mass stars 
present, their lifetimes may well be 1-2 orders of magnitude 
longer. 

For co mpleteness, the Gal a ctic SFR is found to be ~ 2 
Mq yr -1 (|Chomiuk fc Povichl . l201ll ; IDavies et all l201ll ). 
suggesting that a few tens to hundreds of G305 com- 
plexes are analogous to the entire star formation rate of 
the Milky Way. Using results from the Wilkinson Mi- 
crowa ve Anisotropy Probe (WMAP), iMurrav fc Rahman! 
(|2010T ) identify some 31 Galactic HII regions with an ionis- 
ing flux greater than that of G305, of which some 18 WMAP 
sources constitute over half the total Galactic ionising flux. 
Just as the IMF is dominated by the more massive stars 
present, the Galactic SFR is probably dominated by the few 
rigorous star- forming regions present. 



4.4 Alternative Star Formation Rate Indicators 
Within G305 

Clearly the value for the SFR derived from the population of 
embedded massive star-forming regions identified will be an 
upper limit, since we have assumed a power-law slope (i.e. 
Salpeter IMF), and by extrapolating the IMF over a small 
sample of massive stars have overestimated the total mass in 
stars. To this must also be added issues of completeness, and 
assuming a timescale of « 10 5 yrs will be unrepresentative 
of the intermediate to low-mass YSOs present. However the 
value should prove a good assumption for the upper limit, 
that can then be compared to other SFR tracers that are in- 
dependent of the resolved massive stellar population within 
G305. 

What is key here is that we are able to resolve the YSO 
population within Galactic HII regions, such as G305, and 
use both the IMF and stellar timescales to derive a SFR. 
This can then be contrasted and calibrated to extragalactic 
emission tracers, such as the total IR luminosity, to deter- 
mine whether Galactic SFRs are consistent with extragalac- 
tic SFR indicators. We next consider the SFR using tracers 
that are independent of the identified YSO population. Ta- 
ble [3] lists the calculated SFRs for G305 using numerous 
tracers, with reference to each approach. 



44.1 



The Relation Between Star Formation Rate and 
Molecular Cloud Mass 



Recent work by jHeiderman et all (|2010h and 
lLada. Lombardi fc Alve sf (|2010l ) on the star formation 
activity of molecular clouds within 0.5 kpc of the Sun, 
suggest that the star f ormat i on rat e scales linearly with the 
molecular cloud mass. lLadal i| 19921 ) showed that active star 
formation is to be found primarily in high volume density 
regions of molecular clouds, with star formation favouring 
very massive, dense cores. The expectation for the star 
formation activity is thus that there is a tight correlation 
with the amount of high extinction material present in 

mol ecular clouds. 

lLada. Lombardi fc Alvesl (|2010l ) illustrate for a sample 
of local molecular clouds, that by comparison of the cumu- 
lative mass to YSO content as a function of extinction (see 
lLada. Lombard! fc Alvesl <|2010h Fig. 3), a marked minimum 
dispersion of the cumulative mass is found at Ay = 7.3 ±1.8 
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mag. The proposition is that above this minimum, the cloud 
mass is directly related to the star formation activity and 
hence the SFR within the clouds. It is shown that this high 
extinction value correspon ds to an equally high volume den - 
sity of n(H 2 ) ~ 10 4 cm" 3 (|Lada. Lombardi fc Alvesl l201Gf ). 
There is also evidence this linear relation holds for extra- 
galactic molecular clouds. A tight correlation between the 
total IR luminosity and the luminosity of the HCN molecule, 
which itself requires high densities (> 10 4 cm -3 ) to be ex- 
cited t o a detectable le vel, has been found for both Galactic 
cores jWu et al.l . [2005f ). and for a samp l e of n ormal spirals 
and starburst galaxies ()Gao Sz Solomonl l2004h . These stud- 
ies suggest that the linear relation holds for dense interstellar 
gas both on a Galactic and extragalactic scale, underlying a 
physical relation that links star formation and galaxy evo- 
lution. 

From this linear correlation, a SFR is derived of the 
form: 



SFR = 4.6 ± 2.6 x 1(T 8 M . 8 [M yr" 1 ] (10) 

where Mo. 8 corresponds to the cloud mass, in solar masses, 
above an extinction threshold of Ak ~ 0.8 mag, which is de- 
rived from the visual extinction of Av = 7.3 ± 1.8 mag men- 
tioned earlier. 

For each IR source identified by Herschel across the 
G305 complex we determine the physical radius at 250 um, 
as this wavelength offers the optimum combination of signal- 
to-noise, and angular resolution. To estimate the physical 
diameter for each source, we firstly deconvolve the source 
size from the Gaussian beam (|Thompson et all . l2004h : 



©Source — ©Obs 



© 



(ii) 



where &obs i s the source size estimated from the FWHM 
of the Gaussian fitting, and ©s eam is the Herschel beam 
size at 250 um. If we then place each source at a distance 
of 4 kpc, and assume spherical geometry and a uniform den- 
sity, we can determine those sources found to be above the 
critical density threshold. From this we are able to deter- 
mi ne the mass of the dense gas w ithin G305, then using 
the lLada. Lombardi fc Alvesl (|2010l ) assumption that dense 



gas is associated with the star formation activity, we find a 
dense gas mass of ~ 3 x 10 5 Mq ; for comparison, the total 
molecular mass trace d by NH3 is found to be m 6 x 10 5 Mq 
jfflndson et all l201dh ■ We note that this mass is an approx- 
imation, since it is unlikely that the density of each source 
is uniform. 

By combining this mass with Equation (] 10p . we obtain 
a dense gas derived SFR of 0.006-0.02 M yr _1 . The re- 
sult is found to be in good agreement with the embedded 
massive star formation rat e derived earlier, and goe s some 
way toward confirming the lLada. Lombardi fc Alvesl (|2010l ) 
assumption by extension to the star formation activity of 
embedded massive star-forming regions. 



4-4-2 The 70\im Emission Star Formation Rate 

In contrast to deriving a value of the SFR from all identi- 
fied embedded massive star-forming regions, we can also ap- 
proach measuring the SFR by considering the total infrared 
flux (TIR) of the giant G305 HII region. On an extra-galactic 



perspective current SFRs are calculated using a tracer of 
UV photon emission from Y SOs, and also spectral synthe- 
sis models (|Kennicuttl . 1 19981 ). In this case, observations of 
HII regions prove ideal measures of current star formation. 
However, some fraction of the total UV emission will be 
obscured by the presence of dust, thus bolometric IR obser- 
vations of dust (i.e. TIR) will provide an excellent means 
to recover the extinguished UV photon emission, with the 
dust absorption highly peaked in the UV and re-emission 
being in a broad spectral range of mid-to-far-IR |Kennicuttl 
1 1998T ) . The conclusion from this is that the TIR will provide 
the best indicator of SFR obscured by the presence of dust. 

Observationally, an advantage would be the use of a 
single-band star formation indicator, such as UV, Ha, 8 um, 
24 urn etc., however each of these has its own complications. 
In the case of UV and optical lines, corrections need to be 
taken into consideration due to large extinction. Whereas 8 
and 24 um emission strongly depends on the local environ- 
ment, since the abundances of small dust grains that con- 
tribute to their emission depend greatly on metallicity and 
the presence of ion i sing r adiation. 

ICalzetti et~aT1 (120071 ) and iDale et all (120051 ) note that 
8 urn emission makes for an inaccurate SFR indicator since 
there is a large degree of variability of emission in galax- 
ies with respects to SED shape and metallicity. This strong 
variability at 8 um is emphasised in Fig. O where variations 
of an order of magnitude exist between the best fit and good 
fit SED models. This variation can be accounted for by the 
disk inclination to the line of sight for the centrally embed- 
ded object, where the observed flux from a pole-on view can 
be 2 - 4 times greater than a more edge-on viewing angle 
I Whitney et~aU |2003l) . 



Calzetti etall |2005l ) find that the SFR calculated from 



24 um emission itself varies strongly from galaxy to galaxy. 
On a local scale, the ratio of the 24 um luminosity to SFR 
is found to be a reasonably accurate tracer, however when 
applied to other systems, such as starbursts and ultralumi- 
nous infrared galaxies (ULIRGs), the ratio is systematically 
higher. This variation in the 24 um is found to be a factor 
of a few with respects to the observed SEDs, and may be 
grounded in the strong dependence on local galactic condi- 
tions; with ionising stars heating dust to different averaged 
temperatures, the 24 um emissio n will be most sensitive to 
this. However, IDale et all (|2005l ) find that 70 um emission 
may be an accurate monochromatic star formation indica- 
tor, since the 70-to-160 um ratio is found to correlate well 
with local SFRs. 

Recent work bv lLawton et all (|2010l ) have determined 
an accurate monochromatic IR band that best approximates 
the obscured SFR in the Large Magellanic (LMC) and Small 
Magellanic Clouds (SMC), through IR aperture photometry 
of 16 LMC and 16 SMC HII regions, using Spitzer IRAC 
(3.6, 4.5, 8 um) and MIPS (24, 70, 160 um) bands. It is found 
from the IR SEDs of each HII region, that the majority peak 
at around 70 um at all radii (10-400pc) from the centrally 
ionising sources, and that the 70 um emission most closely 
traces the size of each HII region as found using the TIR. 
The conclusion from this is that the 70 um emission is the 
most likely suitable IR band to utilise as a monochromatic 
SFR indicator. 

It has been argued by iKennicutt] (| 19981 ) that the TIR 
is the best obscured SFR indicator available for starburst 
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galaxies. However, dust obscured star formation in HII re- 
gions are found to behave similarly, in that their environ- 
ments are b oth very dusty an d are sites of recent star for- 
mation. The iKennicuttl l| 19981 ) obscured SFR equation is of 
the form: 



SFR = 4.5 x 10" 



[M Q yr" 



(12) 



where Ltib. is the TIR luminosity in erg s , and the value 
4.5 x 10~ 44 is a constant derived from synthesis models, 
with assumptions on the IMF and star formation timescales 
jKennicuttll 19981 ). 

The TIR luminosity in Equation ()12|) can be substi- 
tuted with the averaged 70um luminosity, normalised by 
the TIR, while also applying an IR band spe cific constant. 
The m onochromatic obscured SFR equation of lLawton et al.l 
(|2010l ) is found to be: 



SFR = 9.7(0.7) x KT 44 L(A) [M Q yr -1 ] 



(13) 



where L(A) is the observed luminosity in erg s~ . 

By employing aperture photometry of the whole G305 
region, we obtain the cumulative 70 um flux, f„(A), and from 
this ar e able to calculate the monochromatic luminosity at 
70 um (|Calzetti et all. l2010bh : 



L(X)=4nd 2 (£) f v (\) [ergs 



(14) 



where d is the distance to the G305 complex, in m. 

Using the value found for the observed luminosity at 
70 urn , with Equation (|13|l , we obtain an obscured SFR 
for G305 of 0.002- 0.005 M yr" 1 A similar approach has 
been sugg ested by iLi et all (20101 ) . who also determine a 
monochromatic SFR indicator at 70 urn , yet calibrate their 
SFR tracer not with the TIR luminosity but rather with 
the comb ined Ha, and 2 4 um luminosity. For completeness, 
using the lLi et all (|2010l ) tracer, we derive a SFR of 0.004- 
0.008 M yr" 1 , which is in appr oximate agreement with 
that using the lLawton et all |2010l ) approach. 

We can directly compare these results to that derived 
from the total Lyman continuum photon rate of G305, 
where we find an SFR of 0.002 - 0.004 M yr" 1 (Hindson 
et al. in prep). We note that both these two independent 
tracers are in excellent agreement, however are found to 
a factor of ^ 2 lower than that derived from the embed- 
ded mas sive star-forming p o pulat ion. A similar result is 
found bv lChomiuk fc Povichl (|201ll ), who find that the SFR 
for M17 estimated from both the Lyman continuum and 
24 urn emission is underestimated by a factor of ^ 2 in com- 
parison to the SFR derived from YSO counting. 

This discrepancy may simply be that there are elemen- 
tary differences in the mea surements between G alactic and 
extragalactic observations. lLawton et ail (|2010h note that 
Equation (|13p extends to HII regions measured at projected 
distances of 52 kpc and 61 kpc for the LMC and SMC respec- 
tively. With G305 some ~ 4 kpc distant, the relation between 
the SFR and luminosity at 70 urn may indeed break down, 
with the effects of individual protostars becoming more im- 
portant, due to the larger spatial resolution. It may also be 
the case that the 'steady-state' assumption breaks down in 
this comparison. Though it is normally true that the life- 



time of the region observed is longer than for the individual 
objects for extragalactic realms, this may not hold for ob- 
servations of Galactic regions that tend to be smaller, and 
with shorter dynamical timescales. 

4-4-3 A Galactic - Extragalactic Comparison 

By comparing the derived SFR from numerous tracers, 
shown in Table [3] what is immediately apparent is the dis- 
parity between the rates derived from the resolved stel- 
lar population and those from extragalactic tracers; there 
is lack of consistency between the two, with extragalac- 
tic tracers tending to underestimate the SFR derived 
from resolved Galactic SFRs. This circumstance between 
the two SFR reg i mes h as been noted by several a uthors 
dHeiderman et all. I20l61: lLada. Lombardifc Alvesl . l2010l : 
IChomiuk fc Povichl . 20 111 ), where there a ppears a distinct 
under estimation for Galactic H II regions. iHeiderman et all 
(|2010l ) make the suggestion that this difference may be ac- 
counted for by the inclusion of diffuse gas, in the stan- 
dard Kennicutt-Schmidt relation, that is below the criti- 
cal density threshold for star formation, as suggested by 
lLada. Lombardi fc Alves] l|2010l ). Since extinction maps are 
not readily available to determine the surface density of gas 
in ext ragalactic studies , CO maps are often employed in- 
stead. IHeiderman et all (|2010l ) find that using CO as a gas 
tracer, for a sample of local molecular clouds, leads to an 
underestimate in mass of > 30 % compared to that obtained 
using extinction maps. The result of this would essential! 
push down the estimated SFR from the IKennicuttl (|l99i 
relation, and may go some way in accounting for the dissimi- 
larity between extragalactic regions and more local, Galactic 



However, recent work by iKrumholz. Dekel fc McKeel 

|2012f ) seems to suggest a unified star formation law, with 
objects ranging from both low mass Solar neighbourhood 
clouds through to sub-mm galaxies in agreement with one 
distinct star formation law. What is advocated in this law is 
that the SFR, within a variety of scales, is simply ~ 1 % of 
the molecular gas mass per local free-fall time. This volumet- 
ric approach suggests a local, universal star formation law 
that is applicable to Galactic and extra galactic observations 
(see IKrumholz. Dekel fc McKeel (|2012l ) Fig. 3), bridging the 
gap between the apparent disparity in the two regimes. This 
law is affected solely by local variations, such as the gas 
condition, with more global Galactic/galaxy-scale proper- 
ties, such as the orbital period, having no impact on the 
SFR in so far as they do not change the local properties of 

star-forming reg ions. 

Conversely, lLada et all (|2012l ) conclude that a univer- 
sal star formation law that is applicable from the Milky 
Way through to near-IR selected (BzK) galaxies, is sim- 
ply directed by the amount of dense molecular gas that 
can accumulate within a star-forming region. In the ma- 
jority of situations, only 10 % of the total mass within a 
molecular cloud is at a su f ficient density, n(H2) ^ 10 4 cm -3 
l|Lada. Lombardi fc Alvesl |2010| ). to actively form stars. 
Clearly there is a disparity between the two propo sed uni- 



versal SFR laws; IKrumholz. Dekel fc McKeel (|2012h favour- 
ing g as surface densities and local free-fall times as crucial, 
while lLada et all (|2012l) advocate gas surface densities and 
the fraction of dense gas as the pivotal factors. It therefore 
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Table 3. Calculated SFR for G305 using multiple SFR tracers. 



SFR Tracer 



SFR 
(M yr" 1 ) 



Embedded Massive Stars 0.01 - 0.02 This Paper 

Dense Gas 0.006-0.02 Lada. Lombardi fc Alves (2010) 



UC HII Regions 

Danks 1 & 2 

70 um Emission 

Lyman Continuum 



^ 0.002-0.004 
0.002-0.005 
0.002-0.005 
0.004-0.008 
0.002-0.004 



Hindson et al. (2012) 



Davics 


ct 


al. 


(2012) 


Lawton 


et 


al. 


(2010) 



Li et al. (2010 ) 
Hindson et al. in prep. 



seems that more work is needed to describe the underlying 
nature of a universal star formation law, if such a law is to 
be found. 

What is apparent, when measuring the Galactic SFR, 
is the need for an accurate means to compare the Milky 
Way to other galaxies, in order to allow us to extend the 
more detailed Galactic analysis to other systems and to test 
the discrepancy between the two regimes. Continued, multi- 
wavelength analysis of Galactic HII regions, now including 
Herschel Hi-GAL, will in part aid with this. Through ex- 
tended study across a wide range of star-forming regions, an 
accurate determination of the IMF, and with that the SFR, 
can be achieved, which scaled up from a more local level to 
a global, Galactic level, will allow for the consideration of 
how these crucial properties va ry as a function of en viron- 
ment across the Milky Way fsee lVeneziani et al.l l|2012h for a 
detailed study of the Herschel Science Demonstration Phase 
fields). This should help in a better understanding of how 
the SFR can accurately be measured on both Galactic, and 
extragalactic scales, and lead to a more unified calibration. 



5 SUMMARY 

We have studied the G305 star-forming complex, using Her- 
schel Hi-GAL far-IR data in search of an embedded massive 
star-forming population. In total, we identify some 16 em- 
bedded massive star-forming regions from their associations 
to radio, maser, and RMS counterparts across the region 
(Fig. [7]). From this sample of known embedded massive stars 
we suggest a far-IR colour selection of log (F70/F500) ^ 1 
and log (F160/F350) ^ 1.6 that can be utilised across similar 
regions within the Herschel Hi-GAL survey, to identify em- 
bedded massive star-forming candidates across the Galactic 
plane. 

With our sample of embedded massive stars, we derive 
the present-day SFR for the complex of 0.01 - 0.02 M Q yr~\ 
which is found to be in good agreement with other well 
known m assive star-forming complexes such as the Carina 
comp lex (|Povich et alll201lf ). and M17 (|Povich fe Whitnevi 
2010;). In comparison to other we ll know n extragalactic SFR 
tracers, based on the iKennicuttl (1 19981 ) relation, there is a 
noted discrepancy between the two regimes, with extragalac- 
tic tracers tending to underemphasise the SFR . We note 
however, the use of the lLada. Lombardi fc Alves! (|2010h re- 
lation for SFR that is in good agreement with the resolved 
SFR derived from our study, suggesting the key to the SFR 
is the total amount, and production of, dense gas within 
molecular clouds. 

We find that the Milky Way SFR (|Chomiuk fc Povichl . 



l201ll ; I Davies et all l201ll ) is comparable to tens to hun- 
dreds of G305 complexes, emphasising the fact that 
the Galactic SFR is most likely dominated by simi- 
lar massive star-formi ng complexes. By com bining Her- 
schel Hi-GAL data (|Molinari et all l2010bl) with cur- 
rent and futur e multi-waveleng t h lar ge datasets such 
as GLIMPSE dBeniamin et all 120031). UKIDSS GPS 
' Lucas et all l2008h. MIPSGAL (ICarev et all |2009|) VVV 
Minniti et all |2010|). CORNISH (fPurcell fc Hoarel . I2OIOI ) , 



and MALT90 (|Foster et all |201lJ ) we will be able to study 
massive star formation across the Galactic plane, and how 
the SFR of numerous complexes vary as a function of en- 
vironment. The result could conceivably give rise to global 
star-formation scaling laws that could then allow us to ex- 
tend our detailed knowledge of our Galaxy to more distant 
extragalactic realms. 
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